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Abstract: The Atmosphere Revitalization 

Recovery and Environmental Monitoring 
(ARREM) project was initiated in September of 
2011 as part of the Advanced Exploration 
Systems (AES) program. Under the ARREM 
project, testing of sub-scale and full-scale 
systems has been combined with multiphysics 
computer simulations for evaluation and 
optimization of subsystem approaches. In 
particular, this paper describes the testing and 
modeling of the water desiccant subsystem of the 
carbon dioxide removal assembly (CDRA). The 
goal is a full system predictive model of CDRA 
to guide system optimization and development. 

Keywords: adsorption, flow, thermal transport, 
validation. 

1. Introduction 

Predictive simulation tools are being 
developed to reduce the hardware testing 
requirements of the ARREM project as part of 
NASA’s AES program 1,2 . Although sub-scale 
testing is required to establish the predictive 
capability of the simulations, the much greater 
cost of extensive full-scale testing can be limited 
to that required for the confirmation of analytical 
design optimization studies. Once predicative 
capability is established, geometric 

reconfiguration of a model is usually 

straightforward. A predictive simulation 
capability provides improved understanding of 
complex processes since process conditions 
(temperature, pressure, concentrations, etc.) may 
be examined anywhere in the sorption column. 
Weaknesses in a prototype design can be readily 

identified and improvements tested via 

simulation. Finally, the predictive simulation 
provides a powerful tool for virtual 
troubleshooting of deployed flight hardware. 

Here, we discuss using the COMSOL 
Multiphysics code 3 to model in detail - and 
predictively - experiments that are similar to the 
desiccant subcomponent of a full CDRA system. 

Adsorption in packed fixed beds of pelletized 
sorbents is presently the primary means of gas 
separation for atmosphere revitalization systems. 
However, structured sorbents are emerging as a 



Figure 1 . CBT test stand. 


new approach to sorbent systems. Structured 
sorbents are produced as monoliths, with an open 
structure for airflow, or by fixing sorbents on an 
inert substrate such as paper-like honeycomb 
structures or expanded metal sheets. An accurate 
assessment of structured sorbents and 
comparison with packed bed designs is desirable; 
experimental results so far show unanticipated 
variation in packed bed breakthrough for 
identical beds held under the same conditions. It 
is suspected that small packing irregularities can 
propagate downstream in large beds and impact 
process efficiency. This indicates a margin of 
error inherent in packed bed fabrication and thus 
a likely superiority of structured sorbents for 
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process efficiency and control. This paper 
discusses fully predictive modeling results using 
COMSOL’s Multiphysics code for a 
geometrically simple fixed bed design. Insights 
learned from this work will be used in future 
modeling of the entire CDRA system. 

For the bulk separation of C0 2 and H 2 0, 
temperature changes due to the heat of 
adsorption are significant, requiring the 
simulation of the heat balance equations through 
the beds and the housing, as well as the 
equations for sorption processes and fluid flow. 
For columns with small tube diameter to pellet 
diameter ratios, as encountered in internally 
heated columns, flow channeling along the 
column wall can have a strong influence on 
overall performance. Here, with over a dozen 
pellets per cylinder diameter, 1 -D models should 
prove accurate enough for predictively driven 
system design. 


2. Cylindrical Breakthrough Test 
2.1 Description 

The CDRA requires a water- saving bulk drying 
stage prior to downstream C0 2 removal. The 
primary goal is to continuously remove at least 
80% and up to 100% of water vapor from a 
process air stream. The Cylindrical 
Breakthrough Test (CBT) was constructed to 
compare sorption kinetics for various sorbent 
and sorbate pairs. The tests consist of flowing a 
constant amount of sorbate and carrier gas 
through a fixed bed containing a regenerated (or 
dried out) sorbent. After some period of time 
(the ‘breakthrough time’), the sorbate is detected 
at the bed exit. A plot of the sorbate 
concentration or partial pressure versus time is 
the breakthrough curve. Axially routed 
thermocouples are used to acquire temperature 
curves inside and outside of the bed as well as 



before and after the packed bed. The three 
thermocouples inside the bed are on axis and are 
located 1 inch inside the inlet and exit and in the 
axial middle of the sorbent-filled bed. Absolute 
and differential pressure is measured at the 
column inlet and across the column, respectively. 

The test apparatus, shown in Fig. 1, was 
designed to have relatively low mass (to reduce 
regeneration time) and good axial symmetry (to 
reduce system complexity). The test bed, packed 
with regenerated sorbent pellets, is insulated to 
minimize system heat loss. Mass flow 
controllers are used to blend N 2 , the carrier gas, 
with the desired partial pressure of C0 2 or H 2 0. 
For C0 2 tests, Sable Systems C0 2 analyzers 
provide inlet and outlet C0 2 partial pressure 
readings, while for H 2 0 tests, a Sable Systems 
Dew Point Generator provides humidity control 
and Edgetech Dewmaster dew point analyzers 
provide inlet and outlet dew point measurements. 
A schematic of the entire test setup is shown in 
Fig. 2. The cylindrical column of sorbent has a 
bed length, B , of 16.51 cm and a radius, R B , of 
1.74 cm. The mass of regenerated sorbent in the 
test bed is measured so as to determine the mean 
porosity of the bed. The A1 6061 housing is 
<i=0.065 inches thick and extends for 6 inches 
upstream and downstream of the sorbent. The 
sorbent is held in place inside the housing using 
spring-loaded plates and fine mesh screens. 
Here, we focus on the CBT adsorption tests, 
where the sorbent starts fully regenerated. 
Experimental data from the CBT will be used to 
validate the simulation process so that 
simulation-driven optimization may be used 
alongside conventional design methods to 
perfect CDRA sub-component design and 
testing. 

2.2 Models 


the sorbent pellets, D , the local permeability 
within each cell is then found 4 : 


However, the correct diameter to use in Eq.(l) is 
a volume- weighted averaged such as the Sauter 
mean, but this is not as well-known as D for all 
sorbents, but it is generally larger than D. The 
factor, K.+-, is adjusted to reproduce the measured 
pressure drop through the test bed for a given 
test (so it is not a free parameter). Eq.(l) was 
developed for typical sandstone packing, while 
some sorbents are very spherical and smooth. 
Further, in 1-D the use of a constant porosity in 
Eq. (1) will over-estimate the center-line 
permeability, so values of Kf less than unity are 
typical. 

A 1-D PDE for the interstitial fluid flow 
velocity, u , is derived based on a simplification 
of the Navier-Stokes and Brinkman equations 
while assuming compressible flow; this is 
combined with PDEs for the loading, q, based on 
the Toth equations, the sorbate concentration, c, 
and the temperature of the gas mixture (77), 
sorbent pellets (T s ), and housing walls ( T c ). Note 
that COMSOL flow modules (not used here) all 
solve for the superficial velocity ( [=£u ). The 

interstitial velocity is found via Eq (2), where P 
is the total gas pressure of the carrier gas and 
sorbate mixture, p g is the density of the gas 
mixture, A is the ratio of sorbent area to volume 
(see below), p g is the mixture gas viscosity, and 
is the mean molecular mass of the sorbate. 
Eq.(2) is essentially a 2 nd order Ergun 5 equation 
that reduces to Darcy’s Law in steady state and 
in the limit of a zero velocity gradient and no 
sorption. The 2 nd order Ergun term, proportional 
to the square of the fluid velocity, is somewhat 
akin to a Forchheimer drag term. The — term 
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1-D models of the CBT were constructed using 
the COMSOL Multiphysics code using Domain 
ODEs and DAEs. Only the sorbent-containing 
part of the bed is modeled. A constant porosity, 
t, was used, such that, together with the known 
constant density of the dry sorbent pellets, p s , the 
proper measured total sorbent mass is recovered. 
Using the measured mean particle diameter of 


compensates for the transfer of momentum from 
the overall flow to the stationary sorbent. In 
COMSOL, the entirety of the RHS of Eq.(2) is 
written as a source term in the General Form 
PDE; the pressure gradient cannot be written as a 
flux term or the combined boundary conditions 
become over-constrained. Thus, the boundary 
conditions applied to Eq. (2) are simply 

7*1 = T” I =0. Eq.(2) assumes the ideal 


gas law for the density of the gas mixture, so 
that: 


Pa = 




(3), 


where the mixture’s specific gas constant is 
given by R s =RJM mix . The mean molecular mass 
of the gas mixture is M mix = M N + R g T g c(M a - 
M n )/P , where M N is the molecular mass of the 
pure carrier gas, generally N 2 . 

The second PDE solves for the total gas 
pressure, P, using the continuity equation and 

Eq.(3): 
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In COMSOL, in the Coefficient Form PDE, the 
1 st term of Eq.(4) is a Damping term, the 2 nd is a 
Conservative Flux term, and the 3 rd is an 
Absorption term, while the RHS is a source term 
accounting for the change in total pressure as 
sorption proceeds. The outlet pressure boundary 
condition is of the Dirichlet type, using the 
measured outlet pressure from a given test, P ou t= 
P in -AP. The inlet uses a mass flux boundary 
condition based on the measured constant 
standard flow rate, F, for a given test: 

where Af= 7iR B 2 is the inlet free flow area and P 0 
and T 0 are the reference pressure and 
temperature, respectively, at which the flowrate 
is defined. Thus, the 1 st term of the RHS in 
Eq.(5) is a reference density. Here, P 0 = 1 atm 
and T 0 = 0°C. Since this is a 1-D ‘plug flow’ 
model, there is no gradient of pressure, velocity, 
concentration, temperature, or porosity in the 
radial direction, although out-of-plane thermal 
effects are taken into account (see below). 

The third and fourth PDEs are coupled and 
solve for sorb 


ate concentration, c, and pellet loading, q , 
respectively. The latter applies the Finear 
Driving Force model 6 . Together these two PDEs 
are referred to as the ‘Mass Balance’ equations. 
The General Form PDEs are: 
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where D x is the local time-dependent axial mass 
dispersion coefficient, q* is the equilibrium 
loading from the Toth isotherms, and k m is the 
constant mass transfer coefficient (see below). 
For the concentration, a zero gradient mass 
concentration Constraint: 
c d Pg 
Pg dx 


dc I 

d A =5 


= 0 ( 8 ), 


is used at the exit, while a molar volume flux 
boundary condition is used at the inlet: 

D /fc I _ f_fyg\ _ c 5M m!> I \ _ 

V^Le=0 Pg dx ljr=p dx f=o/ 

EU in C [n (9), 


where u in = is the inlet 

interstitial velocity and c in = P V apin(t)/Rg/T g in(0 5 is 
the sorbate concentration at the upstream inlet to 
the bed. P vap in(0 and T gin (/ L ) are the measured 
time-dependent sorbate partial pressure and gas 
temperature, respectively, for each experiment. 
Since T gin (0 is measured far upstream of the start 
of the sorbent in the bed, the COMSOF inlet 
value is the measured value minus an offset due 
to losses or gains along the uninsulated piping. 
That is, this ensures that T gin (r=0) actually equals 
the T g (£=0,v=0) inside the bed for any given 
experiment. There are no explicit spatial 
boundary conditions required for Eq.(7) since 
there is no loading outside of the sorbent bed. In 
COMSOF, it is necessary to make all of the 
gradient terms of Eq.(6) a flux term, in order to 
apply the inlet flux boundary condition of Eq. (9) 
properly. The above transport PDE, Eq.(6), 
represents Fickian diffusion 6 in conservative 
form and assumes that all mechanical dispersion 
effects are lumped together with molecular 
diffusion in the axial dispersion term. Note that 
together Eqs. (2), (4), and (6) model the transport 
of a concentrated species. That is, M mix is not 
assumed to be constant and the density gradient 
of the carrier gas can have an effect on the 
sorbate concentration. 
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The fifth, sixth, and seventh PDEs are coupled 
and solve for the sorbent temperature, the gas 
temperature, and the wall housing temperature, 
respectively. Together these are referred to as 
the ‘Thermal Balance’ equations. Since the 
sorbent pellets and sorbate plus carrier gas are 
not in thermal equilibrium, the COMSOL Heat 
Transfer module is insufficient and explicit 
PDEs are needed. The General Form PDEs are 
given in Eqtns (10-12). 

All of the terms on the RHS of Eqs.(10-12) 
are source terms in COMSOL, with the 
perimeter terms corresponding to out-of-plane 
convective flux terms. Since u s and u w are 
identically zero everywhere (thus, u=u g ), zero 
flux boundary conditions are used for Eqs.(10 
and 12): 
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In Eqs.(10-12), subscripts s, g, A, and c refer 
to properties of the sorbent, gas mixture, ambient 
environment, and can housing, respectively. Pj = 
2 tzR b , P 0 =2n(R B +d), and A c = n((R B +d) 2 -R B 2 ) are 
the can inner perimeter, outer perimeter, and 
cross sectional area, respectively. Heat capacity, 
thermal conductivity, and thermal transfer 
coefficients are given by c p , k , and h , 
respectively. The heat of adsorption for a given 
sorbate/sorbent pair is given by and is a 
(poorly known) function of q and T s (see Table 1 
below). The ratio of sorbent area to volume 
assumes spherical pellets and is given by A = (1- 
e)6A/B. The outlet boundary condition for 
Eq.(l 1) is given by a zero gradient Constraint: 

ini =°(i3), 

OX l*=F 

while the inlet is given by a flux boundary 
condition: 
dT | 

Tf I — P& ^5 in ( 0 

PX ■jf;= u 


The void fraction is left explicitly in Eq. (14) to 
illustrate that in COMSOL, one must be careful 
to define boundary conditions that are consistent 
with any Conservative Flux terms. In the limit 
of T s =T g , Eqs. (10-14) are the same equations as 
given in a combined Heat Transfer in Solids and 
Porous Media COMSOL module that is set up 
with heat sources and out-of-plane convective 
heat flux nodes. 

The initial conditions for the c and q PDEs 
correspond to equilibrium loading for a small 
initial sorbate partial pressure, typically 1-5 Pa. 
The initial velocity is everywhere, while the 
initial pressure is set to P in -AP *x/B. All 

temperatures are initially linear so as to start with 
the measured (but offset) T fl!n (t = 0) at x=0 and 
V-.ut (t = 0) at x=B. 

The boundary conditions described above for 
P, T, and c are all the measured values from a 
given test, with T and c being time-dependent. 
This is necessary since, as will be seen below, 
the variation from test to test (or even within a 
single test), even using the same packed bed and 
nominal test conditions for flow rate and inlet 
sorbate partial pressure, is often larger than the 
uncertainty in the model. 


2.2.1 Dimensionless Numbers 

Many of the physical parameters used in 
Eqs.(10-12) are determined from correlations 
based on dimensionless quantities. The 
empirical relationships used here are appropriate 
for 1-D models in the regime of the CBT. The 
Schmidt number, Sc, pellet Reynolds number, 
Re, Peclet number 7 , Pe, Prandtl number, Pr, gas 
to sorbent Nusselt number, Nu GS s , and gas to can 
Nusselt number, Nu G c, are given by: 
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In Eq.( 17 ), is 2 for axial ( Pe x ) and 10 for 
radial (,Pe r ) parameters. In Eq.(15), D AB is the 
binary mass diffusion coefficient 10 for either 
water vapor or carbon dioxide in N 2 or air (the 
carrier gas) as a function of P and T g 


2.2.2 Toth Isotherms 

The loading equation uses the Toth isotherm 
relationships: 
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where P vap = c R g T g is the sorbate partial 
pressure and a 0 , b 0 , t 0 , c 0 , and E are Toth 
coefficients for a given sorbent/sorbate pair. The 
coefficients used 12 are listed in Table 1, together 
with the heat of adsorption 12 ’ 13 for each pair. It is 
known that dH depends on loading, decreasing 
significantly in magnitude as equilibrium loading 
is approached as well as when the sorbent is 
colder; the quantitative dependence, however, is 
not well known for many sorbents. A crude 2 nd 
order in loading and 1 st order in temperature fit 
to data for most sorbent/sorbate pairs was done 
at MSFC; only C0 2 on SG is missing. The 
other dH values listed in Table 1 are for a sorbent 
temperature of 12.5°C and a specific loading of 1 


mol/kg. For the silica gel sorbents, Grace Grade 
40 and Sylobead B125, the same Toth 
parameters are used, even though they are 
experimentally derived from the former. The 
same is true for the 5 A zeolite sorbents, Grace 
Grade 522 and RK-38. There are no CBT 
experiments as yet with H 2 0/5A or C0 2 /SG 
systems. 


2.2.3 Material Properties 
Some material properties are from COMSOL’s 
material libraries and some use the correlations 
in the previous section. Temperatures in the 
following expressions are in Kelvin and the 
derived values are in MKS units. 

The axial mass dispersion coefficient used in 
Eq.(6) for the transport of the sorbate through the 
bed is given by: 


uD 

D x = —{ 25). 

P e x 

while for the carrier gas, COMSOL expressions 
for N 2 are used. For the can properties, 
COMSOF functions for A1 6061 are used. Also 
from the COMSOF material library, the heat 
capacity, viscosity, and thermal conductivity of 
water vapor or carbon dioxide, as appropriate, 
were used for the sorbate. 

k gx = CpgPgVx (26) 


The effective axial gas thermal conductivity, 
kgj,, comes from a similarity assumption 6 , which 
may not be valid in all regimes. The resulting 
values are typically as much as 50% larger than 
when using other, more complicated 
expressions 14 ; the impact of different kg X 
expressions on the breakthrough curves 
discussed here is minor. Note that k gx depends 
on the gas mixture heat capacity (see below). 

For the sorbent properties, since they are not 
well known, constant values are used 12,15,16 . 


Table 1 . Adsorption Parameters for Sorbent/Sorbate Pairs 


Sorbate/Sorbent 

System 

mol kg' 1 kPa" 1 

kPa" 1 

E 

K 


K 

m 

kJ mol" 1 

C0 2 /5A 

9.875xl0" 7 

6.761xl0" 8 

5.625xl0 3 

2.700x10'' 

-2.002x10' 

-39.9 

H 2 0/5A 

1.106xl0' 8 

4.714xlO' 10 

9.955xl0 3 

3.548x10' 

-5.114x10' 

-60.5 

C0 2 /13X 

6.509x1 0' 3 

4.884xl0' 4 

2.991xl0 3 

7. 487x1 O' 2 

3.805x10' 

-36.8 

H 2 0/13X 

3.634x10'” 

2.408x1 O' 7 

6.852xl0 3 

3.974x10' 

-4.199 

-56.7 

co 2 /sg 

7.678x10'” 

5. 164x1 O' 7 

2.330xl0 3 

-3.053x10' 

2.386xl0 2 

-40.0 

h 2 o/sg 

1.767xl0 2 

2.787x1 O' 5 

1.093xl0 3 

-1.190xl0' 3 

2.213x10' 

-36.8 


They are listed for the soi 2. Properties of Sorbents as Packed in the CB 


m is less coupled to the flow. 


Sorbent 

Ps 

kgm' 3 

(@26°C) 
J kg 1 K 1 

K 

W m' 1 K' 1 

mm 

mass 

g 

£ 

5A (Grace Grade 522) 

1190 

750 

0.152 

2.22 

125.0 

0.331 

5A (RK38) 

1370 

650 

0.144 

2.10 

119.3 

0.445 

13X (Grace Grade 544) 

1260 

800 

0.147 

2.19 

107.4 

0.457 

SG (Grace Grade 40) 

1240 

870 

0.165 

2.90 

111.7 

0.415 

SG (Sylobead B 125) 

1240 

870 

0.151 

2.25 

127.0 

0.348 


in Table 2. For Sylobead B152, the sorbent 
density and heat capacity are assumed to be the 
same as for Grace Grade 40. Also listed in Table 
2 are the mean pellet diameter, total dry sorbent 
mass, and resulting mean porosity for the sorbent 
packings used in the CBT experiments. 

The gas mixture quantities are calculated 
using a mass-fraction weighting between the 
carrier gas and the sorbate gas. The sorbate mass 
fraction, / a is given by c M a / while the carrier 
mass fraction,/^, is l-/ a . 

The heat capacity of the pellets can change 
significantly as they get loaded. Thus, a mass- 
weighted loading-dependent effective heat 
capacity is used for c ps : 

( C tJ i'-rf 2 £?.££? 2 1 0.CO2 ) + C VSD Ps ) , r - 

v = ~~ — ft; " \ ~ 

yqw {H29C02') -F Ps} 

Eq.(39) typically reduces the pellet heat capacity 
by -10% when the sorbate is H 2 0, but is 
insignificant for C0 2 . 

The thermal transfer coefficients for the gas 
are determined from Nu : 


h sg = Nu es kg/D (2S) 




Nu GC k g 
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(.29), 


while for the ambient-to-can heat transfer, a 
nominal constant small value of h Ac = 0.1 W/m 2 /K 
was used, since the insulation was not firmly 
bonded with the cylinder. As a result, the value 
used for 7^ , the ambient temperature, is not very 
critical, so a constant 19°C value, representative 
of the CBT laboratory, was used in this work. 
However, when available T A was set to the 
measured value for each test, resulting in slightly 
improved late time temperature profiles (see 
below). Also, h E g is typically -150 W/m 2 /K, so 
the gas and sorbent are moderately well 
thermally coupled, while is typically -50 


As a result, in the results discussed below, 
typically, I T g - T s I «1°C and \T g -T c \ »1°C 

3. COMSOL Results 

The only free parameter required to fit the 
CBT data is fc m , the mass transfer coefficient 
used in the LDF model (see Eq.(7)). It is not 
expected to be sensitive to different test 
conditions, such as flow rates, vapor pressure, or 
temperatures. Ideally, it is only a function of the 
sorbent/sorbate pairing, but in practice it can 
vary due to, for example, geometry differences 
when R b is only a few times D, making arbitrary 
predictive applications with the LDF model 
problematic. Here, the same value of is used 
for all CBT experiments of a given 
sorbent/sorbate pairing. Note that drives the 
slope of the vapor pressure rise curve, so for the 
model to be self-consistent, it should reproduce 
that slope; that is, even k tTL is not a true free 
parameter even for a single test of a given 
sorbent/sorbate pairing. The intent is to derive 
fe m for a given sorbent/sorbate pair and use it in 
future predictive modeling of CDRA-related 
systems. It is to be emphasized that once k m is 
determined for a given test of a sorbate/sorbent 
system, the other tests of that sytem (at different 
flow rates and partial pressures) are predictively 
modeled. 

The COMSOL results for the predicted and 
experimentally measured exit temperatures and 
vapor pressures for some illustrative 
sorbent/sorbate pairs are shown in figures 3-6. 
The derived values of used for the various 
sorbent/sorbate systems are listed in Table 3. It 
can be seen that k m varies by a factor of 5, with 
0.002 s" 1 being a typical value. 


Table 3. COMSOL Mass Transfer Coefficients 


3.5 C0 2 on 5A Zeolite RK38 


Sorbent/Sorbate System 

K (s' 1 ) 

Sylobead/H 2 0 

0.002 

Grade 40/H 2 O 

0.00125 

Grade 544/H 2 0 

0.0007 

RK38/C0 2 

0.003 

Grade 522/C0 2 

0.0035 


3.1. H 2 0 on Silica Gel Grade 40 


Two flow rates, 16 SLPM and 8 SLPM, and 
two inlet partial pressures, corresponding to dew 
points of 0.5 °C and 10°C, were run with silica 
gel Grace Grade 40 sorbent. The two flow rates 
correspond to Reynolds numbers of -150 and 
-70, respectively. Figs. 3 and 4 show the 
temperatures and 
vapor pressures 
for the 16 SLPM 
and 0.5°C dew 
point case. The 
COMSOL peaks 
are too low and 
they fall off too 
slowly, 

particularly close 
to the inlet. Also, 
the late time 
temperatures do 
not separate as 
they do in the data. 

The COMSOL 
vapor pressure rise 
matches the data 
well, only missing 
some of the 
complexity of the 
curve shape. This 
is likely due to a 
combination of the 
simplicity of the 
LDF model and 
inaccuracies in the 
dew point data. 

The partial pressure results show that the data at 
late times do not approach the inlet values; the 
reason for this consistent deficiency is unknown 
but may be related to calibration issues with the 
dew point sensors, since it does not occur in the 
C0 2 tests (see below). 


Two flow rates, 16 and 8 SLPM, and two 
inlet partial pressures, 2.5 and 5.0 Torr, were run 
for C0 2 on RK38. Temperature and partial 
pressure comparisons between the experiments 
and the COMSOL models for the 16 SLPM 5 
Torr case are shown in Figs. 5 and 6, 
respectively. As for the SG results, the model 
peak temperatures are too low and fall off too 
slowly. This suggests that issues with the 
thermal transport in general (such as the 
correlation-derived thermal transfer coefficients 
given above) are more likely the problem than 
sorbent- specific properties such as the heat of 
adsorption. Fig. 6 shows that the COMSOL 
model matches the breakthrough curve quite 
well. 


The adsorption capacity, the moles of sorbate 
per unit mass of sorbent that the system can 
adsorb, was measured to be 1.18+0.05 mol/kg. 
In the models, the loading throughout the bed is 
within a few percent of the equilibrium loading 
by the end of the tests. Thus, the calculated 



Figure 3. Temperatures for SG Grade 40 with 16 SLPM and 0.5°C dew point. 
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Figure 4. Water vapor pressures for SG Grade 40 with 16 SLPM and 0.5°C dew point. 



Figure 5. C0 2 temperatures for 5A RK38 at 16 SLPM with 5 Torr inlet 

partial pressure. 


due to too low of an 
activation temperature. 
This would make them 
inaccurate at low C0 2 
loading levels. 

7. Conclusions 


Using COMSOL, 
we have been able to 
derive a predictive 
model for adsorption 
physics in a variety of 
flow rate and sorbate 
partial pressure 

regimes and for a 
variety of sorbents. 
The limitations are in 
general due to the 
experimental data, not 
the model, due to 
variations and 

uncertainties in the 
dew point 

measurements, the 
ambient conditions, and 
unsteady flow conditions. 


However, the model consistently underpredicts 
the temperature peaks and their falloff. This 
could be due to either inaccuracies in the 
isotherms at low loading or the correlations used 
for the heat transfer coefficients. Also, there are 
still deficiencies in the model, such as the LDF, 
which uses a single constant parameter to 
determine the loading rate. Further, the 1-D 
nature of the model is not capturing some of the 
physics due to channeling of the carrier gas near 
the walls of the sorbent bed. Future work will 
extend the PDEs used here to 2-D axisymmetry 
and will include the impact of using a radially 
dependent porosity due to packing. Since no 
purge gas is perfectly dry, another extension of 
this work is to include binary Toth relationships, 
particularly for the C0 2 tests, in order to properly 
capture the competition between H 2 0 and C0 2 . 
The values of fe m found in this work will be used 
in broad applications of sorbent/sorbate systems 
to estimate performance without the need of 
testing. 

8. References 

! Perry, J. L., et, al., "Integrated Atmosphere 
Resource Recovery and Environmental 
Monitoring Technology Demonstration for Deep 
Space Exploration," International Conference on 
Environmental Systems. 

AIAA, San Diego, 2012. 

2 Knox, J.C., et. al., 

"Development of Carbon 
Dioxide Removal Systems 
for Advanced Exploration 
Systems," International 
Conference on 

Environmental Systems. 

AIAA, San Diego, 2012. 

3 comsol, comsol 

Multiphysics®, 2009. 

4 Rumpf, H., and Gupte, A.R., 

"The influence of porosity and 
grain size distribution on the 
permeability equation of 
porous flow", Chemie Ing. 

Techn. (Weinheim), v. 43, no. 

6, p 367-375, 1975. 

5 Ergun, S., Chem. Eng. Prog., 

48, 89, 1952. 

6 D. M. Ruthven, Principles of 
Adsorption and Adsorption 
Processes, Wiley Interscience: 

New York, 1984. 


7 M.F. Edwards and J.F. Richardson, Chem. Eng. Sci. 
24, 607 (1969). 

8 Wakao, N. and Kaguei, S.1982 Heat And Mass 
Transfer In Packed Beds, Gordon and Breach, New 
York. 

9 Li, C.H., and Finlayson, B.A. 1977 Heat transfer in 
packed beds - a reevaluation Chem. Eng. Sci.32 1055- 
1066. 

10 E.L. Cussler, Diffusion: Mass Transfer in Fluid 
Systems, 2nd Edition (Cambridge University. Press, 
London, 1997. 

11 Wang, Y. and Levan, M.D. “Adsorption Equilibrium 
of Carbon Dioxide and Water Vapor on Zeolites 5 A 
and 13X and Silica Gel: Pure Components”, Journal of 
Chemical and Engineering Data, 54(10), 2839-2844, 
2009. 

12 Ritter, J.A., and Ebner, A.D., “Design of an 
Adsorption-Based Carbon Dioxide, Humidity and 
Trace Contaminant Removal System”, Final Report 
for Cooperative Agreement NNM05AA10A, 2008. 
13 ”Davison Molecular Sieves Adsorption Equilibria”, 
Davison Chemical Report. 

14 Yagi, S. and Kunii, D., “Studies on Heat Transfer 
Near Wall Surfaces in Packed Beds”, AIChE J. 6 543- 
546, 1960. 

15 Kay, R., and Pancho, D., “Evaluation of Alternative 
Desiccants and Adsorbents for the 
Desiccant/ Adsorbent Bed”, Honeywell Report 12- 
77742, CAGE 70210, 2013. 

16 Finn, J.E., Ho, E.C., and Knox, J.C., “Progress 
Report on the 4BMS Adsorption Characterization 
Study”, NASA Ames Research Center, 1995. 



Figure 6. C0 2 partial pressures for 5A RK38 at 16 SLPM with 5 Torr 
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